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An  optical  instrument  for  use  in  determining  the  two-dimensional 
vrave  slope  spectrum  of  ocean  capillary  waves  is  described.  The  instru- 
ment measures  up  to  a  35°  tip  angle  of  the  surface  normal  by  measuring 
the  position  of  a  refracted  laser  beam  directed  vertically  upward  through 
the  water  surface.  High  quality  optical  lenses,  a  continuous  two- 
dimensional  Schottky  barrier  photodiode,  and  a  pair  of  analog  divider.- 
render  the  signals  independent  of  water  height  and  insensitive  to  laser 
beam  intensity  fluctuations.  Calibration  is  performed  entirely  in  the 
laboratory.  Sample  records  and  wave  slope  spectra  are  shown  for  one- 
dimensional  wave  tank  tests.  Sample  records  of  two-dimensional  ocean 
tests  are  presented.  A  mechanical  wave  follower  mechanism  was  used  to 
adjust  the  instrument's  position  in  the  presence  of  large  ocean  swell 
and  tides.  Errors  due  to  use  of  a  finite-sized  laser  beam  are  discussed. 
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CHAPTER  I 
INTRODUCTION 


Statistical  descriptions  of  typical  shapes  of  the  sea  surface  and 
its  changes  with  time  are  useful  in  many  studies.  The  sutdy  of  the 
effective  forces  governing  the  interactions  between  wind  and  water  led 
Cox  to  examine  the  statistics  of  small  capillary  waves.   As  Cox  points 
out,  the  shape  parameter  used  in  the  study  of  small,  high  frequency 
waves  is  the  wave  slope.  This  is  so  because  small  wavelength  waves  have 
small  amplitudes  which  become  difficult  to  measure  accurately  as  the 
wavelength  decreases.  However,  these  same  waves  may  have  wave  slopes 
of  considerable  magnitude. 

The  surface  slope  statistics  have  been  useful  in  studies  relating 
to  the  behavior  of  surface  waves  in  the  presence  of  internal  waves,  the 
resonant  interaction  among  capillary  waves,  backseat tering  and  radiative 
transfer  of  electromagnetic  radiation  in  connection  with  remote  sensing 
of  the  state  of  the  sea,  and  relations  between  a  roughened  sea  surface 
and  multipath  effects  in  communications  between  satellites  and  aircraft 

.r.-,  ■       ,,        2-8 
flying  over  the  ocean. 

Optical  methods  of  measuring  waveslopes  seem  to  be  popular  because 
the  .ra.er  surface  is  left  undisturbed  by  the  measurement,  and  the  instru- 
mentation required  to  make  the  measurement  is  basically  simple,  at  least 
in  concept.  Some  of  the  optical  methods  used  in  waveslope  studies  are 
based  on  the  principle  of  measuring  reflected  light.  These  schemes 
include  the  measurement  of  light  reflected  from  an  artificial  source 


and  light  reflected  from  the  sky,  and  the  observation  of  sun-glitter 

0  ?  10  11 
patterns  on  the  open  sea  surface.  '~'  .  '    Other  optical  methods  of 

wave  slope  measurement  are  based  on  the  principle  of  refraction  of  light 
at  a  surface  of  the  water  where  a  discontinuity  in  the  index  of  refrac- 
tion exists.  Measurements  based  on  the  refraction  of  light  have  used 

both  continuous,  extended  light  sources  and  highly  localized  light 

1  1?  13 
sources  in  the  form  of  a  laser  bea.v. .      " 

The  instrument  described  in  this  paper  is  a  second  generation  devel- 

12 
opment  of  the  one  described  by  Tober  et  al.    The  new  instrument  was 

designed  for  and  used  in  both  laboratory  and  oceangoing  field  tests  to  mea- 
sure the  two-dimensional  wave  slope  at  a  point  as  a  function  of  time. 


CHAPTER  II 
DESCRIPTION  OF  THE  INSTRUMENT 


Theory  of  Operation 

The  orientation  of  the  surface  normal  of  an  elemental  area  on  the 
surface  of  a  water-air  interface  may  be  established  by  measuring  the 
angle  of  refraction  of  a  light  ray  passing  through  the  interface .  The 
relationship  between  the  incident  ray  and  the  refracted  ray  is  given  by 
Snell's  Law: 


(l)         n.  sin  ■©:  =  n  sin  ■©- 
i      1    r     r 


where  n,  is  the  index  of  refraction  for  the  medium  containing  the  inci- 
dent ray,  n  is  the  index  of  refraction  for  the  medium  containing  the 
refracted  ray,  and-O:  and -O-  are  the  angles  of  incidence  and  refraction, 
respectively.  By  defining  a  deflection  angle,  <j)  as 

(2)  <j>  =  ^  .. -^ 

and  using  Snell's  Law  as  expressed  in  equation  (l),  the  deflection  angle, 
(P,  may  be  expressed  in  terms  of  n  and-©ras 

(3)  <{)  =  arcsin  (n  sin  ■&-.  )  -  -9-. 


where  n  is  the  relative  index  of  refraction,  defined  as  n  =  n./n  . 

r  r 

Using  the  direction  of  the  incident  ray  as  a  reference  direction,  the 
deflection  angle,  $,  may  he  -used  to  measure  the  refraction  of  the  inci- 
dent ray.  By  aligning  the  incident  ray  with  the  surface  normal  under 
calm-water  conditions,  the  angle  of  incidence,  -fh  ,  may  be  related  to  the 
angular  displacement  of  the  surface  normal  from  the  calm-water  surface 
normal  orientation,  ■©-,  "by 


00      \  -  +  • 


Figure  1  shows  the  geometry  corresponding  to  the  situation  described 

above.  Thus,  substituting  equation  (4)  into  equation  (3).  the  expres- 

.       12 

sion  for  <p   becomes 

(5)  $  =  arcsin  (n  sin  «©■)  -  -d- 

The  inverse  relation  which  expresses  *  as  a  function  of  f   is  readily 
found  to  be 


(6)  -0-  =  arctan 


[sin  $      "J 
n  -  cos  <j) 


Thus  it  is  seen  that  the  measurement  of  the  wave  slope,  •©-,  results  from 
a  measurement  of  the  deflection  angle,  $,   provided  the  relative  index  of 
refraction,  n,  is  known.  Throughout  this  paper  the  value  for  n  is  taken 
as  n  -  1.333. 


direction  of  incident  ray 
and  calm-water  surface 
normal 


plane  of  incidence 


incident  ray- 


Figure  1:  Geometry  associated  with  the  refraction  of  a  ray  of  light  at 
a  water-air  interface.  The  wave  slope  is  defined  as  the 
angle  •©". 


Optical  Design 


The  initial  Gaussian  or  paraxial  design  of  the  optical  components 
for  the  wave  slope  instrument  was  completed  in  late  May,  19?4.  The 
system  was  designed  to  measure  wave  slopes  of-  magnitude  up  to  45  using 
a  United  Detector  Technology  SC/50  Photodiode  as  a  two-dimensional  light- 
to-electrical  signal  transducer.  Matrix  methods  of  ray- tracing  were 

Ik 
employed  to  give  the  first-order  system  design  shown  in  Figure  2.    The 

optical  receiver  portion  of  the  instrument  consisted  of  three  lenses  and 
the  photodiode  transducer. 

In  the  ideal  case,  the  laser  beam  was  assumed  to  be  collimated 
perfectly  and  focused  to  a  point  at  the  focal  plane  by  the  objective 
lens  .  The  position  of  each  ray  in  the  back  focal  plane  of  the  objec- 
tive lens  is  determined  using  a  linear  transformation 


(7)    /  =  *  + 

where  P  is  the  radial  displacement  of  the  focused  laser  beam  from  the 
optical  axis,  a  is  a  constant  which  is  determined  from  the  lens  para- 
meters such  as  surface  curvatures,  index  of  refraction  of  the  lens  mate- 
rial, and  the  thickness  of  the  lens,  and  <j)  is  the  angle  of  deflection  of 
the  laser  beam  from  the  optical  axis.  An  expression  for  (f  in  terms  of 
the  wave  slope,  ■9-,  was  given  earlier  in  equation  (5)«  The  azimuthal 
angle,  "r,  is  unaltered  in  the  transformation.  Thus  it  is  seen  that  any 
light  ray  entering  the  objective  lens  having  a  deflection  angle  (p  and 
azimuthal  direction T  is  directed  to  a  point  (pit)  in  the  back  focal 
plane  of  the  objective  lens.  This  occurs  for  every  ray  entering  the 
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Figure  2:  Paraxial  design  of  the  optical  system. 

The  laser  beam  enters  the  objective  lens  and  is  focused 
130  mm  from  the  back  principal  plane  (BBP).  The  front 
principal  plane  (FPP)  of  the  field  lens  is  located  at  the 
back  focal  plane  of  the  objective  lens.  The  field  lens 
confines  the  laser  beam  to  the  optical  tube  dimensions. 
The  imaging  lens  images  the  field  lens  FFP  onto  the 
photodiode  surface. 


objective  lens,  regardless  of  the  point  of  incidence  on  the  first  sur- 
face of  the  objective  lens.  Since  a  unique  point  in  the  focal  plane 
(fi'f)   corresponds  to  a  unique  ray  orientation  (^YO,  and  since  <f   is 
uniquely  related  to  the  wave  slope,  •£-,  via  equation  (5),  it  is  clear 
that  measurement  of  the  position  the  laser  beam  in  the  focal  plane  (fir) 
gives  unambiguous  information  about  the  orientation  of  the  water  surface 
normal .  Thus  the  problem  of  wave  slope  measurement  is  essentially 

solved  in  the  back  focal  plane  of  the  objective  lens. 

Although  the  necessary  wave  slope  information  was  displayed  in  the 

back  focal  plane  of  the  objective  lens,  it  was  necessary  to  reduce  the 
size  of  the  information-laden  portion  of  the  focal  plane  to  match  the 
dimensions  of  available  photodiode  transducers.  This  was  accomplished 
by  using  a  field  lens  to  confine  the  laser  beam  within  the  dimensions 
of  the  optical  tube  and  an  imaging  lens  which  formed  an  image  of  the 
back  focal  plane  of  the  objective  lens  on  the  active  surface  of  the 
photodiode.  Each  of  the  lenses  shown  in  Figure  2  was  to  consist  of  one 
or  more  lens  elements  to  correct  for  optical  aberrations. 

The  system  described  above  was  never  fabricated  because  a  test- 
schedule  time  constraint  did  not  allow  enough  time  to  have  the  elements 
ground, tested,  and  assembled.  Therefore,  the  design  was  modified  so 
that  off-the-shelf  optical  components  could  be  used.  The  imaging-lens 
requirement  was  filled  by  selecting  a  high  quality  35mm  camera  lens. 
An  objective  lens  was  selected  after  six  weeks  of  testing  a  variety  of 
different  commercially  available  lenses.  The  one  selected  was  a  mili- 
tary surplus  aerial  camera  lens,  a  12",  f/2.5  Aero-Ektar.   It  was  found 
to  be  capable  of  measuring  35  wave  slopes,  displayed  negligible  aberra- 
tion errors,  and  had  a  large  enough  clear  aperture  to  allow  measurements 


to  be  made  with  the  objective  lens  at  a  reasonable  height  above  the 
water  surface.  A  suitable  field  lens  was  located,  but  an  excessively 
long  lead  time  for  delivery  led  to  the  use  of  a  diffusion  screen  in 
place  of  a  field  lens.  The  diffusion  screen  had  the  advantages  of  being 
lighter  in  weight  and  easier  to  mount  in  the  optical  tube,  but  had  the 
disadvantage  of  scattering  most  of  the  light  outside  of  the  imaging  lens. 
The  light  flux  onto  the  photodiode  was  thus  lowered  by  approximately  two 
orders  of  magnitude,  and  additional  electronic  amplifier  stages  were 
necessary  to  produce  a  satisfactory  signal.  The  following  section 
describes  the  components  used  to  construct  the  instrument. 


Description  of  Instrument  Components  ' 

Figure  3a  shows  a  phantom  view  of  the  optical  receiver  portion  of 
the  wave  slope  instrument  corresponding  to  the  original  design  concept. 
Figure  3b  shows  a  view  of  the  final  configuration,  and  Figure  3c  is  a 
photograph  of  the  assembled  unit.  Brief  descriptions  of  each  of  the 
major  components  and  comments  on  their  operating  characteristics  are 
given  in  the  following  paragraphs.  Unless  otherwise  stated,  the  physical 
location  of  the  individual  components  is  indicated  in  Figure  3. 

Objective  Lens 

For  a  given  height,  h,  above  the  water  surface,  the  maximum  measur- 
able deflection  angle  of  the  laser  beam,  (^ax  is  dictated  by  the  clear 
aperture  of  the  objective  lens,  d,  according  to  the  relation 


(8)    <J>max  =  arctan  (l/2(h/d))   . 


Figure  J:     The  optical  receiver  portion  of  the  waveslope  instrument. 

(a)  original  concept,  (b)  final  configuration,  (c)  assembled 
unit.   Components:  laser  beam,  A;  objective  lens,  Bj 
diffusion  screen,  Cj  scattered  laser  light,  Dj  imaging  lens, 
E;  two-dimensional  photodiode,  F;  tube  end  cap,  G;  imaging 
module  mount,  H;  optical  tube,  I;  diffusion  screen  adjust- 
ment screws,  J;  N_  purge  lines,  K;  main  support  collar,  L; 
plate  glass  window,  M;  field  lens,  N. 


A  number  of  commercially  available  lenses  were  tested  on  an  optical  bench 

to  determine  (p  as  a  function  of  the  ratio  (h/d).  A  military  surplus 

aerial  camera  lens,  an  Aero-Ektar,  f/2.5,  12"  focal  length  lens,  was 

selected  as  being  the  most  suited  to  the  intended  purpose  of  the  wave 

slope  instrument .  A  special  lens  element  support  was  fabricated  for  the 

lens  to  reduce  the  total  weight  from  18.75  lbs  to  11  lbs.  Figure  4 

shows  the  maximum  measurable  deflection  angle,  &       .  as  a  function  of 

»max 

the  objective  lens  height  above  the  water  surface,  h.  Scales  indicating 
the  maximum  measurable  wave  slope, -0-  ,  and  the  ratio  (h/D),  where  D 
is  the  maximum  diameter  of  the  lens  barrel  at  the. plane  of  the  first 
element  of  the  objective  lens,  are  also  shown  in  Figure  4. 

The  shape  of  the  back  focal  surface  of  the  objective  lens  was  deter- 
mined by  a  series  of  Foucault  knife-edge  tests.  The  surface  showed  con- 
siderable curvature  at  large  displacements  from  the  optical  axis,  but  a 
plot  of  the  radial  displacement  of  the  focal  point  of  the  laser  beam  as 
a  function  of  the  deflection  angle,  <f),   was  linear  to  a  good  approxima- 
tion. The  results  indicated  that  the  linear  transformation  expressed  in 
equation  (r/)  was  accurate  to  within  0.3%  or  less  and  a  Gaussian  or  par- 
axial model  of  the  objective  lens  was  indeed  a  good  approximation.  This 
notion  was  reinforced  later  by  the  results  of  the  calibration  tests  and 
the  measurement  of  errors  due  to  optical  aberrations. 

Diffusion  Screen 

The  diffusion  screen  was  used  in  place  of  a  field  lens  because  a 
satisfactory  lens  could  not  be  procured  at  the  time  of  instrument 
fabrication.  The  diffusion  screen  served  to  scatter  light  from  the 
focused  laser  beam,  thus  allowing  a  sufficient  quantity  of  light  to 
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Performance  curve  for  the  objective  lens.  Symbols  are 
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enter  the  imaging  lens  module  and  activate  the  photodiode.  A  number  of 
different  types  of  translucent  and  frosted  plastic  sheets  were  teste- 
find  a  satisfactory  screen  material.  Scattering  lobes  for  eight  differ- 
ent potential  screen  materials  were  measured  and  evaluated.  Best  results 
were  obtained  using  a  one-eighth-inch-thick  sheet  of  clear  acrylic  which 
was  roughened  to  a  uniformly  frosted  appearance  by  grinding  both  surfaces 
with  a  ^-00-grit  abrasive. 

The  diffusion  screen  was  placed  at  the  focal  plane  of  the  objective 
lens  using  a  set  of  adjustable  mounting  tabs.  The  mount  was  made  adjust- 
able rather  than  rigidly  fixed  to  allow  for  accurate  and  optimal  adjust- 
ment of  the  position  of  the  screen  in  case  the  instrument  calibration 
tests  indicated  that  such  an  adjustment  was  necessary.  Subsequent 
calibration  tests  showed  that  the  initial  placement  of  the  screen  was 
satisfactory  and  no  further  adjustment  was  necessary. 

Imaging  Lens  Module  -  Optics 

The  imaging  lens  module  consisted  of  a  35mm  camera  lens,  the  two- 
dimensional  Schottky  barrier  photodiode,  and  the  electronic  amplifier 
and  signal  conditioning  components.  These  items  were  mounted  in  a  common 
aluminum  housing  and  the  entire  module  somewhat  resembled  a  35mm  camera 
in  appearance.  Figure  5  shows  front  and  rear  photographs  of  the  imaging 
lens  module. 

An  interference  filter  was  attached  to  the  camera  lens  to  reduce  the 
effect  of  background  illumination  on  the  output  signal  from  the  photo- 
diode. This  filter  passed  the  He-Ne  laser  line  at  6328  Angstroms  and 
had  a  bandwidth  at  the  half  power  points  of  100  Angstroms.  A  100X  band- 
width was  selected  because  the  laser  light  scattered  from  the  edges  of 
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Figure  5  s  The  imaging  lens  module;  front  view,  (a),  and  hack  view,  (h), 


the  diffusion  screen  impinged  on  the  surface  of  the  interference  filter 
at  relatively  large  angles  of  incidence  arid  would  not  have  been  passed 
by  a  much  narrower  bandpass  filter.15  The  imaging  lens  was  a  Super 
Takumar  multi-coated,  55mm  focal  length,  f/i.4  camera  lens. 

Imaging  Lens  Module  -  Electronics 

The  laser  beam  position  sensing  transducer  was  a  United  Detector 
Technology  continuous  Schottky  barrier  Photodiode  Model  SC/50  with  a 
3.56  cm  x  3,56  cm  active  area,  shown  in  Figure  6.  Similar  devices  have 
been  used  as  radiation  tracking  transducers,  and  descriptions  of  mathe- 
matical models  for  the  lateral  photoeffect  are  available  in  the  litera- 
ture.     The  photodiode  may  be  operated  in  either  a  photoconductive 
mode,  where  a  reverse  bias  voltage  is  applied  to  the  diode  junction,  or 
a  photovoltaic  mode,  where  no  reverse  bias  is  applied  to  the  diode  junc- 
tion. In  the  photoconductive  mode  the  depth  of  the  depletion  region 
associated  with  the  p-n  junction  is  physically  greater  than  that  occuring 
in  the  photovoltaic  mode  of  operation.  The  photoconductive  mode  features 
a  lower  junction  capacitance  and  higher  junction  resistance  which  gives 
rise  to  a  much  higher  frequency  response  capability,  but  also  it  gives 
rise  to  a  greater  absolute  value  of  noise  due  to  the  reversed-bias  shot 
noise.  The  high  value  of  the  junction  resistance,  approximately  68 
megohms  at  11  volts  reverse  bias,  allows  greater  input  light  intensities 
to  strike  the  diode  surface  with  no  damage  due  to  excessive  junction 
currents.  The  photovoltaic  mode  features  a  higher  junction  capacitance, 
lower  junction  resistance,  and  a  reduced  noise  level  due  to  the  absence 
of  the  reversed-bias  shot  noise.  The  frequencies  dealt  with  in  the 
present  application  were  low  enough  to  make  the  effect  of  the  junction 


Ji* 


JiWla^^^Wl^^ 


= 


i"»'"l.  i-'.i.  i.iJTT.I.i  .l.i.l.i.i-T  „„i,.i.  J"..,  u\u\''^ 


Figure  6  :  Two-dimensional  Schottky  Barrier  Photodiode, 


capacitance  negligible.  Both  modes  of  operation  were  tested,  and  the 
photovoltaic  mode  appeared  to  be  more  stable  during  calibration  tests. 
Therefore,  all  of  the  subsequent  tests  were  conducted  using  the  photo- 
diode  in  a  photovoltaic  mode  of  operation. 

Figure  ?  shows  a  view  of  the  imaging  lens  module  with  the  hinged 
door  opened  to  expose  the  electronic  amplifiers,  the  regulated  power 
supply  and  the  analog  dividers.  Originally,  commercially  available 
amplifiers  were  used  to  amplify  the  photodiode  output,  but  serious  " 
signal  drift  problems  required  re-trimming  of  the  amplifier  stages  every 
15-30  minutes.  This  made  practical  application  of  the  instrument  impos- 
sible .  The  amplifier  shown  in  Figure  7  was  a  first-generation  improve- 
ment assembled  using  low-drift  amplifier  components.  This  amplifier 
allowed  operation  of  the  instrument  for  periods  of  up  to  six  hours  with- 
out incurring  any  serious  errors  due  to  electronic  amplifier  drift. 
Components  for  a  second-generation  improved  amplifier  were  ordered  but 
not  delivered  in  time  to  be  used  in  any  of  the  tests  conducted  and 
reported  on  in  this  paper. 

Figure  8  shows  a  block  diagram  of  the  electronic  amplifier  and 
related  signal  conditioning  components.  The  photodiode  produces  four 
current  signals  through  the  four  ohmic  contacts  to  the  n-type  lattice 
material .  Four  transconductance  amplifiers  convert  the  current  signals 
to  amplified  voltage  signals.  The  sum  and  difference  amplifiers  increase 
the  gain  of  the  signals .  The  analog  dividers  normalize  each  difference 
signal  to  the  sum  signal;  the  quotients  are  the  x-  and  y-axis  outputs 
which  are  insensitive  to  input  light  intensity  fluctuations  at  the 
photodiode .  The  circuits  for  the  x-  and  y-  channels  are  virtually  iden- 
tical, and  are  shown  in  the  schematic  diagram  in  Figure  9. 
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Figure  7  :  Electronic  amplifiers. 
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Figure  8  :  Block  diagram  of  the  electronic  amplifier. 
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Figure  9  :  Schematic  diagram  of  one  amplifier  channel.  Approximate 
values  for  resistors,  R,  and  capacitors,  c,  are:  R.  ,  R.  = 
3.5  Kft,  Hi  =  .6  Mil,  E2  =  .5  Ha  ,  R3  -  10  Kn,  R4  -  3.5  Kn.  , 
B5=5U,  R6-  IKft,  R7  =  50  Ko-  ,  R8  -  50  Kit  ,  c  =  1000  af, 


Low  drift  analog  dividers,  manufactured  by  Analog  Devices,  were 
used  to  perform  the  signal  normalization  with  respect  to  the  input  light 
intensity.  The  dividers  were  trimmed  with  external  potentiometers  in 
the  manner  suggested  by  the  manufacturer.  It  was  found  that  normaliza- 
tion accuracy  was  within  1%  when  the  denominator  voltage  was  in  the 
range  -8  v  to  -4  v.  This  performance  was  somewhat  degraded  as  the  denom- 
inator voltage  increased  to  -10  volts  and  also  an  error  of  Jfo  was  ob- 
served as  the  denominator  dropped  to  -2  v.  The  normalization  error 
increased  very  rapidly  as  the  denominator  voltage  fell  below  -1.5  v. 
Therefore,  a  useful  operating  range  was  taken  to  be  any  laser  intensity 
that  would  produce  an  analog  divider  denominator  voltage  in  the  range 
-2  v  to  -10  v. 

The  analog  dividers  served  as  a  very  convenient  means  of  compensa- 
ting for  variations  in  the  laser  light  intensity.  This  was  especially 
useful  during  calibration  of  the  instrument ,  since  the  analog  dividers 
served  as  real  time  computers  performing  the  signal  normalization.  The 
normalization  could  have  been  performed  more  accurately  and  over  a 
greater  range  of  laser  beam  light  intensity  fluctuations  by  recording 
digitized  voltages  of  the  sum-  and  difference-amplifiers  and  performing 
the  division  on  a  digital  computer.  However,  this  method  would- have 
been  more  cumbersome,  more  time  consuming  and  would  have  caused  a  serious 
delay  in  the  development  of  the  instrument.  The  effectiveness  of  the 
analog  dividers  in  compensating  for  laser  beam  intensity  fluctuations 
is  illustrated  by  the  oscilloscope  traces  in  Figure  10.  Figure  10a 
shows  60  Hz  fluctuations  in  the  sum-amplifier  outputs  due  to  powering 
the  laser  from  a  rectified  60-cycle  power  supply.  The  x-axis  denomi- 
nator is  shown  as  the  upper  trace  and  the  y-axis  denominator  is  shown 
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Figure  10:  Illustration  of  the  effectiveness  of  the  analog  dividers 

to  produce  a  normalized  signal  voltage.  Typical  fluctuations 
in  an  un-normalized  signal,  (a);  the  normalized  signal,  (b). 


as  the  lower  trace .  The  numerator  signals  exhibited  the  same  type  of 
fluctuations  with  roughly  the  same  amplitudes  as  those  shown  for  the 
denominators.  Figure  10b  shows  the  result  of  normalizing  the  numerator 
signals  to  the  denominator  signals  with  the  analog  dividers.  Notice  the 
60  Hz  fluctuations  are  virtually  eliminated. 

The  entire  amplifier  was  found  to  be  very  sensitive  to  small  changes 
in  the  supply  voltage  to  the  amplifier  components.  A  regulated  +12  v 
power  supply  was  incorporated  to  maintain  a.  constant  supply  voltage  to 
within  5  mv.  at  room  temperature,  75   F,  and  to  within  +  20  mv  at  tempera- 
tures ranging  between  45  F  and  95  F. 

Optical  Tube 

The  optical  tube  was  fabricated  using  0.025"  aluminum  sheet,  type 
2024-T.3.  The  sheet  was  rolled  and  spot  welded,  and  then  sealed  with 
room-temperature-vulcanizing  silicone  rubber.  The  external  surface  was 
coated  with  a  titanium  dioxide  epoxy  resin  to  provide  extra  protection 
from  the  elements  and  also  to  reflect  sunlight  and  aid  in  maintaining 
a  lower  internal  tube  temperature.  The  internal  surface  of  the  tube  was 
coated  with  a  flat  black  paint  to  minimize  reflections  of  scattered 
light.  The  main  support  collar,  shown  in  Figure  3i  was  an  annular  ring 
machined  from  acrylic  plastic  stock.  The  support  collar  was  designed  to 
carry  the  entire  weight  of  the  optical  receiver  without  loading  the  thin 
aluminum  tube  in  a  lateral  direction.  This  minimized  the  chance  of  tube 
flexure  and  errors  introduced  by  distortion  in  the  optical  alignment  of 
the  internal  components. 


Oceangoing:  Support  Structure 

The  oceangoing  support  structure  In  its  final  configuration  is 
shown  in  Figure  11 .  The  structure  was  designed  to  withstand  loading  by 
six-foot  ocean  waves  with  deflections  that  would  produce  optical  errors 
no  greater  than  O.jfc.     Laboratory  tests  at  three  to  four  times  the 
design  loading  showed  that  deflection  errors  were  on  the  order  of  +  C.jfc 
to  +  0.7%  with  resonant  frequencies  for  various  deflection  modes  ranging 
between  1.8  and  3«8  Hz.  This  indicated  that  the  flexural  errors  encoun- 
tered during  ocean  operation  would  be  less  than  the  0.3$  design  value. 

The  laser  carriage  had  two  sets  of  adjustment  screws  to  aid  in 
establishing  the  optical  axis.  The  optical  receiver  carriage  featured 
a  double-gimbal  support  ring  to  carry  the  weight  of  the  optical  tube. 
This  made  optical  alignment  in  the  field  a  relatively  simple  task, 
usually  requiring  no  more  than  15-3°  minutes.  A  bulls-eye  type  of  bubble 
level  was  mounted  on  the  upper  horizontal  member  of  the  support  struc- 
ture to  aid  in  adjusting  the  optical  axis  to  the  vertical.  The  level 
was  calibrated  to  provide  a  measure  of  the  optical  axis  deviations  from 
the  vertical  to  +  3  in  any  direction. 

The  entire  oceangoing  support  structure  was  designed  to  be  mounted 
on  the  servo-driven  boom  of  a  wave-follower  mechanism.  The  wave-follower 
mechanism  was  designed  and  built  by  the  Department  of  Coastal  Engineering 
at  the  University  of  Florida.  It  is  a  semi-portable  structure  that  may 
be  mounted  on  the  ocean  floor  in  thirty  feet  of  water  as  shown  schemati- 
cally in  Figure  12.  The  hydraulically  driven  boom  has  a  dynamic  range  of 
+  2  meters  about  the  mean  water  level  and  operates  such  that  an  attached 
instrument  package  maintains  a  relatively  constant  height  above  the 
water  surface. 
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Figure  11 i  The  oceangoing  support  structure. 
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Figure  12  :  Sketch  of  the  wave  slope  instrument  during  ocean  operation. 


CHAPTER  III 

THEORETICAL  PERFORMANCE  AND 

OPERATING  CHARACTERISTICS  OF  THE  INSTRUMENT 


Introductory  Remarks 

Calibration  of  the  instrument  consists  of  measuring  the  amplif ied 
output  voltage  signals  for  each  of  the  two  orthogonal  axes  of  the  two- 
dimensional  photodiode  as  functions  of  the  laser  beam  deflection  angle, 
(j),  and  the  azimuthal  angle ,  ^  The  deflection  angle  ,  fy,    is  then  related 
to  the  wave  slope,  -0?  via  equation  (6)  and  the  orientation  of  the  surface 
normal,  as  shown  in  Figure  i,  may  be  determined.  If  the  signal  normali- 
zation process  discussed  earlier  is  to  be  performed  on  a  digital  computer, 
then  four  voltages  must  be  recorded  for  each  value  of  $   and  t  ,   i.e.  the 
sum-amplifier  output  and  the  difference-amplifier  output  for  each  of  the 
two  orthogonal  axes  is  required.  If  the  analog  divider  is  used  to  per- 
form the  normalization,  then  only  two  voltages,  the  analog  divider 
outputs,  are  required  for  each  value  of  <f   and  }*  .   Calibration  using  the 
analog  dividers  was  by  far  the  more  convenient  method  of  calibration 
during  the  development  of  the  instrument,  but  it  was  subject  to  uncer- 
tainties in  calibration  due  to  some  sensitivity  of  the  analog  dividers 
to  the  value  of  the  denominator  signals. 
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Calibration  Using  the  Analog  Dividers 

Figure  13  shows  a  photograph  of  the  calibration  set  up  in  the 
optics  laboratory.  A  2  raw  polarized  HeNe  laser  was  mounted  at  right 
angles  to  the  optical  bench  which  was  used  to  support  the  optical 
receiver  portion  of  the  wave  slope  instrument.  The  laser "beam  was 
directed  through  a  Polaroid  filter  to  a  right  angle  prism  used  to  reflect 
the  beam  parallel  to  the  axis  of  the  optical  bench.  The  axes  of  the 
optical  bench  and  the  reflected  laser  beam  were  aligned  prior  to  mount- 
ing the  instrument  in  an  optical  bench  cradle.  The  .cradle  was  then 
adjusted  such  that  the  optical  axis  of  the  instrument  coincided  with  the 
path  of  the  laser  beam.  The  Polaroid  filter  could  be  rotated  to  vary  the 
light  intensity  into  the  instrument  and  thus  simulate  attenuation  of  the 
laser  beam  due  to  absorption.  The  prism  was  used  to  simulate  refraction 
at  water  surface.   Changes  in  the  deflection  angle,  (>,  could  be  simulated 
"by  rotating  the  prism.  This  was  accomplished  by  mounting  the  prism  on  a 
leveled  turntable  that  could  "be  rotated  about  a  vertical  axis.  The 
accuracy  of  the  turntable  movement  was  approximately  0.02°.  Changes  in 
heights  above  the  water  surface  could  be  simulated  by  sliding  the  instru- 
ment cradle  along  the  optical  bench. 

Figure  14  shows  a  typical  plot  of  the  calibration  data.  The  x-axis 
and  y-axis  analog  divider  output  voltages  are  denoted  by  0  and  0  , 
respectively.  Solid  lines  connect  data  points  corresponding  to  f  =  con- 
stant and  broken  lines  connect  data  points  corresponding  to  (j)  =  constant. 
The  points  were  generated  by  first  setting  the  null  signal  position,  and 
then  stepping  the  simulated  deflection  angle,  $,  in  2°-increments  across 
the  entire  surface  of  the  objective  lens.  The  output  voltages  from  the 
analog  dividers  were  read  using  a  digital  voltmeter  and  recorded.  The 
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Fig\ire  13  :  The  calibration  set  up  in  the  optical  laboratory. 
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Figure  lh\   A  typical  plot  of  calibration  data  obtained  using  the 
analog  divider.  Q  and  Q    have  units  of  volts. 


instrument  was  then  rotated  approximately  15  about  the  optical  axis  and 
the  procedure  repeated. 

A  major  concern  in  the  application  of  the  instrument  is  the  repeat- 
ability of  calibration  data.   It  was  found  that  at  a  constant  temperature 
the  calibration  data  were  consistent  to  within  2-3%-  However,  a  change 
in  temperature  of  10-15  led  to  calibration  variations  on  the  order  of 
10-20%.  Nearly  all  of  this  variation  could  be  accounted  for  by  temper- 
ature dependent  changes  in  the  offset  voltage.  The  amplified  offset 
voltages  as  a  function  of  temperature  are  shown  in  Figure  15.  The 
temperatures  were  measured  inside  the  optical  tube  at  the  physical 
location  of  the  amplifiers  with  the  optical  tube  sealed  in  its  ocean- 
going configuration.  A  gas-tube  thermometer  with  a  dial  indicator  was 
used  to  make  the  temperature  measurements.  The  thermometer  and  its 
location  relative  to  the  amplifiers  is  shown  in  Figure  l6.   Calibration 
data  for  various  optical  tube  temperatures  were  taken  in  the  laboratory. 
In  the  field,  the  tube  temperature  during  the  test  was  monitored  so  that 
an  appropriate  set  of  calibration  data  could  be  applied  to  the  resulting 
test  data.  After  correcting  for  temperature  changes,  the  maximum  un- 
certainty in  the  calibration  was  estimated  to  be  on  the  order  of  ±1%. 


Evaluation  of  Optical  Aberrations 

Optical  aberrations  result  when  rays  from  a  common  point  in  the 
object  space  enter  the  lens  at  different  locations  and  fail  to  come  to 
a  common  point  of  focus  at  the  location  of  the  image  point  determined 
from  Gaussian  or  paraxial  approximations.  Any  rays  which  deviate  from  the 
Gaussian  approximation  are  thought  of  as  some  sort  of  aberration. 
Lengthy  treatments  of  aberrations  may  be  found  in  any  of  several  good 
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Figure  15  :  Offset  voltages  measured  with  respect  to  the  amplifier 

outputs.  Numerator  and  denominator  voltages  are  denoted 
by  n  and  d,  respectively.  Subscripts  refer  to  axes  of  the 
photodiode.  Error  bars  indicate  the  scatter  of  data 
points;  the  broken  line  segments  represent  temperatures 
at  which  no  data  points  were  generated. 


Figure  16  :  Photographs  showing  the  location  of  the  optical -tube 
thermometer  relative  to  the  amplifiers  in  the  imaging 
lens  module. 


reference  books  on  optics. 

The  effects  of  optical  aberrations  e-ssociated  with  the  objective 
lens  were  evaluated  by  an  indirect  method.  Since  optical  aberrations 
may  affect  both  the  location  of  a  point  in  the  image  plane  and  the 
sharpness  of  the  image,  and  since  any  significant  aberrations  would 
lead  to  measurable  variations  in  the  instrument's  output  signals,  the 
output  of  the  instrument  itself  was  used  to  estimate  the  severity  of 
optical  aberrations. 

The  laser  beam  was  treated  as  a  "fat  light  ray."  This  test  ray  v;as 
allowed  to  enter  various  portions  of  the  objective  lens  while  keeping 
the  deflection  angle,  §,   and  the  azimuthal  angle,  j,  constant.  During 
these  variations  the  output  signals  were  monitored  to  detect  any  signifi- 
cant changes.  In, practice,  the  angle  r  was. held  constant  by  avoiding 
any  rotation  of  the  optical  tube  after  the  optical  axes  of  the  instru- 
ment and  the  laser  beam  were  aligned  on  the  optical  bench  as  shown 
previously  in  Figure  13 .  The  deflection  angle,  j),  was  then  varied  such 
that  the  laser  beam  moved  over  the  entire  surface  of  the  objective  lens. 
The  x-  and  y-axls  output  voltages  were  recorded  and  plotted  as  a  function 
of  <{>.  The  distance  between  the  prism  and  the  objective  lens  was  varied 
keeping  T  constant  and  the  deflection  angle,  f>,  was  again  varied  through 
its  full  useful  range.  Again  the  x-  and  y-axis  voltages  were  recorded 
and  plotted  as  a  function  of  ®.  The  entire  procedure  was  repeated  a 

third  time  at  still  a  different  value  for  the  separation  between  the 
prism  and  the  objective  lens  of  the  instrument.  This  process  allowed 

the  beam  with  a  specific  orientation,  (|>,r)  to  enter  the  objective  at 

different  radial  locations  on  the  surface  of  the  lens.   If  aberrations 

had  been  important,  the  data  points  corresponding  to  common  values  of 


$   and  T  would  have  changed  noticeably.  This  was  not  the  case,  however. 
The  voltages  corresponding  to  common  values  of  j   and  j   were  repeatable 
to  within  the  value  of  the  noise  in  the  analog  divider  outputs.  Thus, 
it  was  concluded  that  the  effects  of  optical  aberrations  were  negligibly 
small  and  completely  masked  by  the  noise  in  the  electronics. 

The  Frequency  Response  of  the  Instrument 

The  frequency  response  of  the  instrument  was  determined  by  physically 
chopping  the  laser  beam  with  a  rotating  mechanical  aperture.  The  laser 
beam  was  positioned  to  produce  a  null  signal  output  from  the  analog  divid- 
ers when  the  beam  was  allowed  to  enter  the  objective  lens.  Blocking  the 
laser  beam  with  a  mechanical  stop  caused  the  denominator  signals  to  the 
analog  dividers  to  go  to  zero,  thus  the  divider  became  saturated  at  the 
maximum  divider  output  voltage.  By  switching  the  laser  on  and  off  at 
high  frequencies  the  frequency  response  of  the  entire  system,  from  photo- 
diode  to  divider  output,  was  evaluated.  The  results  showed  that  the 
frequency  response  of  the  entire  system  was  flat  to  frequencies  well 
above  ^+00  Hz.  Since  all  of  the  data  gathered  with  the  instrument  were 
filtered  to  eliminate  frequencies  in  excess  of  100  Hz  before  digitizing 
the  analog  signals,  it  was  concluded  that  the  frequency  response  of  the 
instrument  was  flat  throughout  its  useful  operating  range. 


wavelength  Response  Characteristics 

Cox  has  shown  that  the  finite  dimensions  of  the  spot  on  the  water 
through  which  all  rays  must  pass  on  their  way  to  the  detector  form  a 


"low-pass"  filter  to  the  signal.  Measurable  wavelengths  must  satisfy 
the  inequality  A  ^6.8r  where  r  is  the  radius  of  the  spot  size,  and  A 
is  the  resolvable  wavelength.  This  resolution  criterion  is  applicable 
to  systems  in  which  the  light  intensity  is  uniform  across  the  spot.  The 
following  analysis  considers  the  case  where  the  light  intensity  varies 
across  the  spot  in  the  form  of  a  two-dimensional  Gaussian  distribution 
function. 

A  first  approximation  to  the  wave  slope  instrument's  response  to  a 
finite-sized  laser  beam  assumes  that  the  photodiode  is  sensitive  only  to 
the  centroid  of  the  laser  beam  focused  on  the  diffusion  screen.  This 
is  a  reasonable  assumption  because  the  spot  on  the  diffusion  screen  is 
usually  small  and  well  defined  for  wavelengths  satisfying  the  Cox  criteria, 

This  assumption,  plus  the  use  of  the  linear  relationship  between  the 
radial  displacement  of  a  ray  on  the  diffusion  screen,  ?  ,  and  the  deflec- 
tion angle,  f,   as  given  in  equation  (?)  allows  the  analysis  to  be  simpli- 
fied. It  follows  that  P  -     $ ,   where  <j>   is  defined  as  the  deflection  angle 
corresponding  to  refraction  at  a  plane  surface  which  results  in  a  centroid 
located  at  F. 

The  variation  between  f>,   the  measured  deflection  angle  and  <fot   the 
deflection  angle  of  the  central  ray  in  the  laser  beam  was  examined  first. 
Next,  the  fact  that  the  photodiode  is  not  strictly  sensitive  to  the  laser 
beam  centroid  was  dealt  with  and  an  error  associated  with  the  spreading 
of  the  light  intensity  distribution  about  the  centroid  was  estimated. 

Errors  Due  To  A  Finite-sized  Laser  Beam 

The  theory  of  opeation  presented  earlier  was  based  on  the  refraction 
of  a  single  ray  of  laser  light.  For  the  single  ray  case,  the  wave  slope, 


©-,  is  related  to  the  deflection  angle,  j,   by  equation  (6): 

sin  $  "} 


(6)    -0-~  arctan 


n  -  cos  $ 


where  n  is  the  relative  index  of  refraction  for  water. 

In  the  case  of  a  laser  he am  of  finite  dimensions,  the  beam  may  he 
considered  as  a  bundle  of  rays  distributed  over  a  finite  region  about 
the  optical  axis.  Here,  the  optical  axis  is  placed  at  the  centroid  of 
the  laser  beam  as  it  emerges  from  the  laser  cavity.  In  general,  rays 
penetrationg  the  water  surface  at  different  locations  will  encounter 
different  wave  slopes  and  hence  give  rise  to- different  deflection  angles. 
A  measure  of  the  wave  slope  at  the  optical  axis  is  still  desired,  but 
it  becomes  necessary  to  work  with  a  mean  deflection  angle,  0,  rather 
than  a  single  deflection  angle  at  the  optical  axis.  A  value  for  the 
measured  wave  slope,  -0-  ,  is  defined  as 


(9)   «"   =  arctan 


m 


o 


Since  </)  will  generally  vary  from  the  value  of  the  deflection  angle  at 
the  optical  axis,  4  ,  the  measured  wave  slope,  -9-  ,  will  generally  be 
different  from  the  true  wave  slope  occurring  at  the  optical  axis,  -&- 
Therefore,  an  uncertainty  or  error  in  wave  slope  measurement  may  be 
defined  as 


(10)     £  =  -e-  -  -e- 

m    o 


It  is  desirable  to  describe  the  behavior  of  the  error,  £,  as  a 
function  of  tne  laser  beam  parameters  and  the  configuration  of  the  water 
surface  at  the  test  site.  This  is  a  problem  of  considerable  complexity, 
but  a  first-order  estimate  of  the  behavior  of  the  error  may  be  made 
following  the  introduction  of  a  few  simplifying  assumptions.  These 
assumptions  are 

(1)  Geometric  optics  are  assumed  to  be  adequate  for  a  satis- 
factory description  of  the  system,  i.e.  diffraction  effects 
are  considered  negligible  and  ray  tracing  techniques  may 

be  employed; 

(2)  The  coherence  properties  of  the  laser  beam  are  assumed  un- 
important, thus  interference  effects  are  negligible  and 

a  light  intensity  distribution  rather  than  an  electromag- 
netic wave  amplitude  distribution  is  associated  with  the 
laser  beam. 

(3)  The  laser  beam  below  the  surface  of  the  water  is  assumed 
to  be  perfectly  collimated,  i.e.  consisting  of  a  bundle 
of  perfectly  parallel  rays. 

(^4-)  The  light  intensity  distribution  in  the  laser  beam  is 
assumed  to  be  a  Gaussian  distribution  centered  about  the 
optical  axis. 

(5)  The  surface  of  the  water  is  assumed  to  be  accurately  de- 
scribed by  a  linear  combination  of  sinusoidal  plane  waves 
of  varying  wave  numbers  and  with  varying  directions  of 
propagation . 


The  fifth  assumption  allows  the  surface  amplitude  to  be  described 
as 
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(11)  z(x,y,t)    =    SNjSinCk.r    -  Wj*  +    £   ) 

where  the  subscript  "i"  refers  to  the  i  -component  of  an  infinite  series 
of  plane  waves,  k  is  a  three-dimensional  wave  number,  r  is  a  radius 
vector  measured  perpendicular  to  the  optical  axis,  w  is  the  radian 
frequency  of  the  wave  component,  t  is  the  time,  and  §  is  an  arbitrary 
phase  offset  measured  at  time  t=0.  The  x-,  y-dependence  is  introduced 
by  the  use  of  a  cylindrical  coordinate  systems 

(12)  z  =  z 

x  =  r  cos  "^ 
y  =  r  sin  "4* 

where  J   is  an  azimuthal  angle  measured  with  respect  to  the  positive 
x-axis.  Figure  17  illustrates  the  configuration  for  the  i  -wave  com- 
ponent propagating  in  the  T   azimuthal  direction.  The  Gaussian  laser 
beam  is  centered  about  the  optical  axis  at  (x=0,  y=0)  and  is  directed 
in  the  positive  z-direction. 

To  simplify  the  analysis,  consider  a  single  wave  component  with 
wave  number  k,  where  the  cartesian  reference  frame  has  been  rotated 
such  that  the  wave  component  travels  along  the  y-axis,  and  %   has  been 
set  to  zero.  Under  these  conditions,  the  surface  is  described  by 

(13)  z  -  a  sin  (ky  -  wt) 
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Figure  17  :  Configuration  for  the  i  -wave  component  and  the 
Gaussian  laser  beam. 
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where  the  subscript  "i"  has  been  dropped.   This  type  of  analysis  is 
consistent  with  that  of  Cox  in  estimating  the  wavelength  response 


.  ..   1 
characteristics . 


The  wave  slope  of  this  simplified  model  is  then 


dz  . 

(14)    -9-=  =  ak  cos  (ky  -  wt ) 

dy 


In  cylindrical  coordinates 

(15)  -e-=  ak  cos  (kr  sinClK)  -  vrt) 
The  deflection  angle,  j>,  is  found  by  Snell's  Law, 

(16)  sin(tf>  +©)  -  n  sin  £-, 

to  be 

(17)  </)  =  arcsin(n  sin(-&))  -  "6" 

which  becomes,  upon  substituting  equation  (15)  into  equation  (l?), 

(18)  <f(r,f)  =  arcsin  \_n   sin(ak  cos(kr  sin(f)  -  wt))J 

-  ak  cos(kr  sinC^")  -  wt)   . 

The  deflection  angle,  tf(r,f),  is  specified  at  each  point  (x,f)   through 
out  the  xy-plane  by  equation  (18).  The  Gaussian"  laser  beam  weights 


each  of  these  points  by  a  factor 


(19)    p(r,t) 


IT 


0" 


where  the  factor  ^  was  introduced  to  make  the  weight  function  normal- 
ized in  the  sense  that 


-too 


(20)       ^p(r^)  dA-  =  1 


The  expectation  value  of  $,  or  the  mean  value  of  j)   is  f     and  is 
defined  by 


OD  Z\T 


(21>     $=    <$>  -    ^|(r,t)  p(r,f)  r  dr  df   . 

o  O 

Substituting  equations  (l8)  and  (19)  into  equation  (21 )  gives  the  ex- 
pression 

~    W    2  2        r 
(22)     $  -  \(  -^-e_<rr  arcsin  1  n  sin(ak  cos(kr  sin(>) 

00  -  wt))lr  dr  dt 

00  ZIF  x   2  2 
-  ((  X  e"  °"r  ak  cos(kr  sin(t)  -  wt)  r  dr  if    . 

o  \ 

Referring  back  to  equation  (17)  and  Snell's  Law,  it  is  seen  that  the 
first  integral  in  equation  (22)  has  no  obvious  physical  significance, 
but  the  second  integral  in  equation  (22)  is  clearly  the  mean  value  of  the 
wave  slope  relative  to  the  Gaussian  weight  function  supplied  by  the 

laser  beam  intensity  distribution.   Integration  of  the  second  integral 

.  .   20 

is  relatively  straightforward,  giving 


<»  iff  «2 

(23)        W  £  ak  e_<rr"  cos(fct  sin  ("*)  -  wt)  r  dr  dT    = 

ak  cos(vt;  e 

Since  ak  cos(wt)  is  the  slope  evaluated  at  the  optical  axis,  the  weight- 
ed mean  slope  is  simply  the  slope  at  the  optical  axis  weighted  hy  the 
factor  exp(-k  /k<r). 

Integration  of  the  first  integral  in  equation  (22)  could  not  be 
accomplished  in  closed  form  so  that  a  numerical  scheme  was  necessary  for 
evaluation . 

The  numerical  integration  scheme  used  an  approximation  to  the  Gauss- 
ian weight  function  as  illustrated  in  Figure  18.  The  central  cylinder 
and  each  cf  the  annular-ring  segments  carried  equal  flux  and  were  thus 
equally  weighted.  These  equally  weighted  segments  were  then  divided  by 
a  number  of  radials  such  that  each  of  the  resulting  sub-segments  were 

equally  weighted.  The  equally  weighted  sub-segments  were  then  replaced 

f 
by  equally  weighted  delta-functions.  The  delta-functions  were  located 

at  the  centers  of  each  sub-segment,  with  the  center  being  defined  by 
a  mean  radial  distance  and  a  mean  azimuthal  angle  for  each  of  the  sub- 
segments.  The  approximation  shown  in  Figure  18  provides  a  l6-point 
approximation  to  the  Gaussian  distribution.   In  practice.  360-point 
and  1800-point  approximations  were  used. 

Note  that  the  solution  to  the  second  integral  of  equation  (22)  as 
expressed  in  equation  (23)  is  a  function  of  three  parameters:  ak,  the 
maximum  wave  slope;  wt,  the  phase  of  the  plane  wave:  and  the  ratio  k/<r  . 
This  suggested  a  change  in  the  radial  integration  variable  for  the  first 
integral  of  equation  (22)  from  r  to  rcr=y  .  After  making  this  variable 
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Figure  18  :  Approximation  to  the  Gaussian  intensity  distribution 
used  in  the  computer  error  analysis. 
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change  and  re-defining  the  locations  of   the  delta-functions  in  the 
Gaussian  weight  function  approximation,  the  mean  deflection  angle, $  , 


became 


(&  n      1  ,.  r 

(24)       J)    ~j2:5:    —7"  S  (j>-   TJ  arcsin   In  sin(ak  cos(k/cr 
(ij 

j?     sin(l^)  -  wt))|f    -  ak  cos(wt)  e" 

Equation  (24)  *as  then  evaluated  on  an  IBM  37°  digital  computer  for 
various  values  of  the  parameters  ak,  wt,  and  k/j-  .  The  same  computer 
program  also  evaluated  -&-   as  expressed  in  equation  (9)  and  the  wave 
slope  error,  €,  as  defined  in  equation  (lO). 

At  this  point  it  is  convenient  to  digress  momentarily  to  clarify 
the  significance  of  <T  and  k/<r   .  The  laser  beam  parameter,  <T  ,  is  a 
measure  of  the  radius  of  the  Gaussian  intensity  distribution  given  in 
equation  (19).  The  radius  of  a  Gaussian  laser  beam  usually  refers  to  the 
l/e2  point  of  the  beam,  i.e.  the  radius  at  which  the  intensity  reaches 
approximately  0.14  of  the  maxium  intensity.  Calling  this  radius  rQ, 
it  is  seen  from  equation  (19)  that  <T   may  be  defined  as 

(25)        <r  =  Y2Vro 

and  has  units  of  reciprocal  length.  Note  that  since  (T  is  a  function  of 
r  ,  it  is  also  a  function  of  distance  from  the  output  window  of  the  laser 
cavity  for  a  diverging  beam.  • 


A  plane  wave  of  wave  number  k  has  an  associated  wavelength  given 
by 

(26)      X-  2TT/k  . 

Thus  it  is  seen  that  the  ratio  of  the  water  wavelength  to  the  laser 
beam  radius  may  be  expressed  as 

2TT  ^.^ 

(2?)     */r_ 


o 


Y21  (k/<r  )       (k/r 


Equation  (27)  is  plotted  in  Figure  19- 

The  parameter  k/<r  was  convenient  to  use  in  the  computer  analysis, 
but  wavelength-to-laser  beam  radius  ratios  are  difficult  to  visualize 
using  k/ff     as  a  parameter.  Therefore,  whenever  k/<r  is  used  as  a 
parameter,  the  conversion  to  units  of  Vr0  is  indicated. 

The  computer  calculated  wave  slope  errors  indicated  that  the 
greatest  absolute  error  occurred  at  phase  wt=0,  i.e.  at  the  point  of 
greatest  wave  slope.  Figure  20  shows  the  error  in  wave  slope  for  various 
values  of  the  parameter  ak  plotted  as  a  function  of  parameter  k/cr  . 
In  all  cases,  the  error  is  associated  with  the  worst-error  case  occurring 

at  wt=0. 

Figure  21  shows  the  percent  error  In  wave  slope  measurement  for 
ak^Jr  5°  and  akjfc:  45°  plotted  as  a  function  of  the  parameter  k/(T  .  Again, 
these  values  correspond  to  the  worst-error  case  of  wt  =  0.   It  is  seen 
that  a.  10fc-   error  in  wave  slope  measurement  occurs  at  k/<T  =0.65 
(A/r0  =  6.33)  for  ak*  5°,  and  \/<r-   O.76  (A/rQ  =  5-84)  for  ak^  hf. 
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Figure  19  :  X/rQ  as  a  function  of  parameter  k/g-  .  Additional  abscissas 

are  scaled  to  show  the  percent  error  and  the  error  measured 
in  degrees  for  the  case  of  a  maximum  wave  slope  of  45°. 
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Figure  20  •  Error  in  wave  slope  measurement , £ ,  as  a  function  of 

parameter  k/<r.  at  wt  =  0,  for  various  values  of  parameter  ak, 
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Figure  21  s  Percent  error  in  wave  slope  measurement  as  a  function  of 
parameter  k/<r,  for  wt  =  0  and  ak  4?  5  ,  k5   . 


Taking  k/g-  =  0.7  (tyrr  =  6.34)  as  a  representative  value  correspond- 
ing to  a  lOfo  error  at  wt  =  0,  the  behavior  of  the  error  as  the  phase  of 
the  wave  changed  from  zero  degrees  to  90  was  examined  for  various  values 
of  the  wave  slope  parameter,  ak.  The  computer  evaluations,  of  the  error, 
"based  on  a  360-point  Gaussian  approximation,  are  plotted  in  Figure  22. 
The  most  noteworthy  features  in  Figure  22  are  the  noticeable  inflection 
points  in  the  curves  for  slopes  >  30  ,  and  the  actual  change  in  sign  of 
the  error  for  values  of  slope  approaching  k-0   .  The  sign  change  was  found 
to  be  a  physically  observable  phenomenon  and  not  an  artifact  introduced 
by  the  computer  analysis.  The  error  in  wave  slope  measurement  was 
found  to  be  satisfactorily  modeled  by  the  expression 


(28)  £(wt)  =  ^(ak,  k/r)  cos  wt  +  £2(ak,  k/j-)  cos  3wt 

where  the  coefficients  €  and  €0  are  functions  of  the  parameters  ak 
and  k/tf* .  For  certain  values  of  these  parameters,  the  coefficient  £ 
may  become  greater  than  £  ,  thus  shifting  a  substantial  portion  of  the 
error  from  a  frequency  w  to  a  frequency  3w.  For  example,  a  steep  sin- 
usoidal wave  having  a  slope  of  ^5  ,  when  measured  with  a  Gaussian 
laser  beam- with  radius  r  "=  0.1 5?  \,   would  produce  a  4.3^  error  at  radian 
frequency  w  and  a  k  .J/o   error  at  a  radian  frequency  3w.  The  result  is  an 
erroneous  indication  that  the  wave  of  radian  frequency  w  contains  a 
superposed  ripple  component  of  radian  frequency  3w.   Growth  of  the  co- 
efficient £2  is  very  rapid  as  the  slope  exceeds  30°,  but  for  slopes  of 
35  or  less,  greater  than  80?2  of  the  error  is  associated  with  the  coeffi- 
cient €.     ,  the  "in  phase"  error  component.  Since  the  instrument  described 
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Figure  22  :  Error  in  wave  slope  measurement,  6  ,  as  a  function  of 
wave  phase  for  various  values  of  maximum  wave  slope 
and  for  k/f  =0.7  (  tyr  -  6.3*0- 


in  this  paper  is  not  suitable  for  measuring  slopes  in  excess  of  35  at 
reasonable  heights  above  the  water  surface,  the  problem  of  errors  shift- 
ing to  higher  frequencies  may  be  considered  of  negligible  importance . 

The  quantity  defined  by  the  symbol  ^-    should  be  taken  as  an  error 
rather  than  a  correction  factor  to  be  applied  to  the  wave  slope  data, 
unless  the  wave  number  spectrum  of  the  wave  field  is  accurately  known. 
This  is  true  because  the  orbital  velocity  or  the  fluid  velocity  associ- 
ated with  the  gravity  waves  destroys  the  one-to-one  correspondence 
between  the  observed  signal  frequency,  f ,  and  the  wave  number,  k, 
as  is  discussed  presently. 

A  simple  wave  field  consisting  of  a  single  sinusoidal  plane  wave 
propagating  through  the  test  site  v^puld  produce  a  wave  slope  signal  of 
a  single  frequency,  f ,  related  to  the  radian  frequency,  w,  by  the 
expression 

(29)    f  =  w/2TT  . 

Phillips  gives  the  radian  frequency  of  a  wave  as 


(30)    w  =  (gk  +  V  k3)l/2 


where  g  is  gravitational  acceleration,  k  is  the  wave  number  of  the  wave, 

and  K  is  a  constant  determined  from  the  surface  tension  and  the 
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density  of  the  water.    He  also  shows  that  the  wave  number  of  a 

short  wave  changes  in  the  presence  of  a  long  wave  according  to 

the  expressions 


(3i)    k  =  k0  ">  *Y. 
(32)  6k  =  ko(l)/c) 

wnere  k  is  !    '  e  number  of  the  short  wave  in  the  presence  of  a 
long  wave,  kQ  is  the  wave  number  of  the  short  wave  in  the  absence  of 
a  long  v:ave,  6  k   is  the  observed  change  in  The  wave  number  of  the  shori 
wave,  U  is  the  fluid  velocity  of  the  long  wave,  and  c  is  the  phase 
velocity  of  the  long  wave.   It  is  readily  found  that  an  approximate 
expression  for  the  radian  fr.   '  ■  -y  of  a  short  wave  in  the  presence 
of  a  long  wave  5 s 

,         1/2 

(33)  w~(gkQ4  Vko3  -  gkou/c+     3lfk0^/c)         . 

The  observed  signal  frequency  may  then  be  written  as 

1/2 

(34)  f  *  (l/2ir)(gko  +  if  V  +  gkoU/c  -  3^ko  u/c) 

At  the  crest  of  a  long  wave,  U  is  a  maximum  positive  value  and  the 
obserV(  __  .,-,    .  f   juency  for  the  short  wave  3s  greater  than  the 
freqlle  lc;        (d  when  no  long  wave  is  present.  When  no  long  wave 
is  present,  U   0  ai  I  equation  (3*0  degenerates  to  equation  (30). 
:  trough  location  of  the  long  wave,  U  is  a  maximum  negative 
•     ,  jrved  frequency  is  less  than  that  occurring  when 
3  Lon<   -  e  is  p] 


Determining  the  fluid  velocity  ,  U,  in  the  presence  of  severa] 
randomly  superposed  long  waves  of  varying  wavelengths  becomes  complex 
and  intractable  for  practical  purposes.  Thus,  the  one-to-one  corres- 
pondence between  an  observed  signal  frequency  and  a  unique  wave  number 
as  expressed  i   i     .  ions  (29 J  and  (30 )  is  destroyed  by  ths  presence 
of  long  waves.  Surface  currents,  whether  generated  by  the  action  of 
wind  on  the  water,  tides,  or  other  mechanisms, also  complicate  the 
problem  by  introducing  additional  fluid  velocity  components. 

The  error  In  the  wave  slo,         icr.t,  as  discussed  above,  is 
not  the  only  source  of  error  due  to  a  finite-sized  laser  beam,  although 
it  -,-,y   be  considered  dominant.   A  secondary  error,  b,    Is  associated 
with  the  type  of  sensor  used,  i.e.  the  continuous,  two-dimensional, 
Schottky  barrier  photodiode.  This  secondary  error  results  from  the 
fact  that      I    iiode  '  >  n<  I   truly  sensitive  to  the  centroid  of 
-the  light  intensity  distribution  striking  the  surface  of  the  photodiode 
That  is,  a  delta-function  distribution  and.  a  Gaussian  distribution, 
each  having  the  same  centroid  location  and  carrying  the  same  amount 
of  flux,  will  generally  produce  different  output  signals  from  the 
photodiode.  Taking  the  delta-function  distribution  as  a  reference 
distribution,  the  error,  &  ,   may  be  defined  as  a  fractional  error 

(35)        £  =  (S  -  S   J/S  -. 

where  2  is  the    ]  it  signal  for  a  general  distribution  of  light  on 
the  surface  cF   the  photodiode  and  S  >f  is  the  output  signal  for  a 


delta-function  distribution,  it  being  iu  Lea   od  bhat  both  distributions 
have  the  same  centroid  location. 

A  conservative  estimate  to  the  magnitud  of  5  for  the  case  of  a 
Gaussian  1  se:   sam  of  radium    interacting  with  a  sinusoidal  plane 
wave  of  wavelength  'X  -  6.3**  rQ  and  a        wave  slope  of  35  was 
found  to  be  6  -  -0.019,  i.e.  a  -1.93  error.  Referring  back  to 
Figure  21,  it  is  seen  that  the  percent  error  due  to  the  error  €   is 
approximately  -1096  for  the  case  A  =  &.;;!j  rQ  and  slope  =  35  .  Thus  . 
the   6 -error  accounts  for  more  thar  :   of  the  total  error.  For  the 
case  of  a      of  35°  and  \   -•>  rQ,  the  $  -error  increases  to 
approxii  tel         id  the   £~erroi  ii        to  approximately  -\%. 
Thus  it  is  seen  that  the  influence  of  the  tf-orror  becomes  more  and  more 
apparent  as  th«   iti<  X/rQ    '  s*«if  cor. tiibuting  significantly  to 
the  total  error  associated  with  the  usu  of   a  finite-sized  laser  beam 
as  the  wavelength,  A,  oache      Limension  of  the  laser  beam. 


■ '■■':■_.    [V 


The  instrument  was  used  in  a  series  of  one -dimensional  wave  tank 
tests  conducted  at  the  University  c^  Florida's  Wave  Tank  Facility  in 
March,  1975-  The  laser  was  mounted  horizontally,  outside  of  the 
tank  roughly  35  cm  below  the  mean  water  level,  The  laser  beam  was 
directed  through  a  pia      is  wind  .;  bo  an  underwater-mounted  prism 
where  the  beam  was  deflected  upward  along  the  vertical.  The  wave  slope 
instrument  was       L  with     ■  ■.,'■  8  cm  above  the  mean  water 

level.  Optical  alignment  war.  performed  under  calm-water  conditions  with 
the  x-  and  y-axis  channels  directed  • .,  the  crcsswir.d  and  downwind  direc- 
tions, re       ily. 

Figure  ?J   sho»rs  samples  of  the  downwii   :h   lei  signal  for  the 
cases!   no       Leal  wave,  wind  speed  of  11.9  m/sec;  2  second,  2.5  cm 
hanical  wove,  wind  speed  of  11.9  v./sec;   and  2  second  5  cm 
amplitude  mechanical  wave  wit},      i  ■:    B.2  m/cec.   In  all  cases 

the  upper  time  series  is  the  downwind  channel  output  signal  and  the 
lower  I       '  es  is  a  capacitance  gauge  ou       Lch  gives  an  indica- 
tion of  the  phase  of  the  wave.   Me       L  waves  and  wind  go  from  .right 
to  left.  Seine  features  to  note  are 

(i)  tl        pn  -nt  of  a  4  Hz  wave  structure  due  to  the  action 

wave, 
(2)   the  h-   H;;  wind        .   .,  ,  .    lergoes        rent  frequency 
■  o  h.'-'.   Hz  at  the  crest  of  bhe  2.5  c~!  mechanical 
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're      wave  tank  tests, 


wave,  ana  decreases  to  j.6   Hz  in  the  trough, 

(3)  the  frequencj  i;:xt.ulation  of  the  wind  generated  nave  over  ■ 
the  5  cr=  mechanical  wave,  with  frequencies  approaching 

10  Hs  at  the  crest  and  '+>  Hz-  in  the  trough, 

(4)  the  concentration  of  ripple  activity  on  the  downwind  face 
of  the  wind  generated  waver,  as  indicated  by  the  80-100  Hz 
components  in  the  wave  slope  instrument's  time  series 
record, 

(5)  maximum  wave  slopes  of  11  -!■-•'  in  the  case  of  the  8.2  m/sec 
wind  speed,  and  20-23  in  the  crr.se  of  the  11.9  m/sec 
wind  si 

Phillips  has  shown  that  the  wave  number  and  the  wave  amplitude  of  short 

21 

waves  increase  at  the  cresl  oi  long    .       de  raase  in -the  trough. 

This  gives  rise  to  greater  wave  slope*,  wii        observable  frequencies 
at  the  crest,  of  the  lent,  wave  and  lesser  wave  slopes  with  lower  fre- 
quencies at  the  trough.  Tl    a  record  for  the  5  cm  mechanical  wave, 
as  shown  in  Figure  23,  is  ei  .irelj  con:       with  Phillips'   prediction. 

At  high  wave  slopes  the  laser. beam  would  deflect  off  the  objective 
lens  and  a  signal  dropout  occurred.   Dropout  was  detected  by  a  zero 
denominator  signal.  Minimal  dropout  occurred  during  these  tests  and 
there  are  no  examples  of  dropout  in  the  sample  records  shown. 

Figure  24  shows  the  frequency  spectra  of  the  .  ean  square  slope  for 
the  three  sample  records  plotted  on  leg-log  scales.   The  energy  units 


ctra  presented  were  prov:c  d  Kr.  Allan  Reece  of  the 
Depart                 Lneering  at  th  v  of  Florida.  The 

spectre  am  computer  program 

develo]  ison,22 


10 


;c~2 


A.. 


i 


.. 


I 
(radians  ) 

ir"" 


10  -' 


io-< 


L0 


J 


0      10 


4  1 


■i  \  -*►    5 


10 


'.::. 


10 


10 


10 


10"^ 


ENERGj 


(radians  -sec) 


•^  -*>  *   10 


-t 


10' 


:  Frequ<        tra  of  the  mean  square  slope, 


2 

for  the  2.5  cm  and  5.0  an  mechanical  wave  cases  are  in  radians  -sec 

and  the  appropriate  ordinate  appears  to  +.h<?  right  of  the  figure.  The 
spectrum  for  the  case  of  no  mechanical  wave  vis  multiplied  by  the  fre- 
quency in  order  to  compare  with  the  spectrum  of  Cox  which  is  also 
shown  in  the  figure.1  The  units  of  these  latter  cases  are  radians  and 
the  appropriate  ordinate  appears  to  the   left  of  the  figure.  Note  the 
increase  in  the  0.5  Hz  component  as  the  mechanical  wave  amplitude  in- 
creases.  Also  note  the  decrease  in  the  <';  Kz  component  as  the  mechanical 
wave  ;       •  increases.        ittLy,  this  is  due  to  a  frequency  spread- 
ing effect  as  discussed  earlier  in  connection  with  the  modification  of 
short  waves  by  the  fluid  velocity  associated  with  long  waves.  The 
reason  for  the  order  ot  magnitude  difference  between  the  Cox  spectrum 
and  the  present  spectrum  is  still  under  active  investigation.  Some 
ible  reasons  for  the  difference  are  the  presence  of  a  light  oil 
slick  in  the  University  of  Florida  Wave  Tank;  the  Cox  fetch  was  2.14  m 
while  the  present  tests  were  conduct  1  at  a  fetch  of  9-1  m;  Cox  was 
sensitive  to  light  intensity  fluctuations  and  also  made  linearizing 
assumptions  concerning  the  light  intensity  distribution  across  the 
snot  and  the  modulation . of  the  intensity  due  to  changing  wave  slopes; 
mm  and  sub-mm  ripples  may  have  affected  the  present  measurements  due 
tc  excessive  spreading  of  the  refracted  laser  beam. 


CHAP]1! 
OCEAN  Tv— 


■ear.  wavp  i  made  on  May-  29?-  19 ^5 »  at  a 

tier,  off  the  coast  of  Marii     ,  Florida.  /  shrimp  boat  carried 
the  auxiliary  power  supply  units  and  the  >       siisition  station.  The 
wave  follower  sei  in  thii        of     r.   A  standpipe 

■  ..  -  '  the  wave-follower  boom  to 

facilitate  the  it..     iti  n  i       .  al  of  instrument  packages.  Figure 
2'  ;'  *aph  of  the  i   Imp  boat's  position  relative  to  the  test 

£J  te  . 

nt  carriea  separate  hi  t   ..  p  >wer  packs  for  the  elec- 
tronic amplifii  ion  Model  80-?s  submersible 
laser.  Wave  clop;'  si    :  rere  hard-wired  to  the  data  acquisition 
center.  Figur   •'■'  shows  the  instrument  in  operatic        the  ocean 
test.  I:   immi    :       test  site  was  protected  from  the  effects  of  sun- 
g] ittcr  fcy  an  opaque  plastic  s 

records  of  1  shown  in  Figure  27. 

(a)  and  (b)        y-  and  x-axis  outputs,  respectively.  Trace 
(c)  is  the  x-.  which  was  monitc  ■■   to  ide: 

■  .  Trace  (d)  is  a  wave  staff  record  indicating 

Ions  of  I      :m  ocean  swell. 

;  for  several  reasons; 

1  i   ary  tests 
a  standpipe  the  instrument  was  snapped  by 

rents  at 


>  1 


Figure  25  :  Photograph  of  's  position  relative  to  the 


cccrji  be: 


Figure  26  .  the  instrument  in  operation  during  the  ocean  test, 
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Figure  27   :   Sample  records  of  ocean  test  data 


(2)  The  test  site  v.  cluttered  wJ  '.,   interference  during  the 

est ,  Swells  ref  ]       :-n  the  shrimp  boat,  and  the 
test  site  was  directly  on  Lhu  lee  side  of  the  shrimp  "boat, 

(3)  The  staudp'  ladder  was  too  close  to  the 
test  site        Ij  frequently  slapped  the  ladder  causing 
serious  waveie       '..icr.z. 

(-; )  Portions  of  the  record  displayed  considerable  dropout  due 
to  fish  in  the  laser  bean.  On  several  occasions  a  group 
of  20  -  y.)   small  fish  toofc  tarns  filing  through  the  laser 
Fish  dro]  identified  by  long  periods 

of  zero  denominator  si, 
(5)  ''■■■  cost  of  com]  it  r  tim       ■  d  to  analyze  the  ocean 
data  was  cons  Sine      Lrect  application  could 

ie  to  items  (l)-  (*l)  above, 
;   ■ . seed. 


SUMMARY 


instrumeni  ibed  here  has  show::  several  oig;;ific:u.;  -Havar.cec: 
in  obtaining  continuous  time  records  of  vivo  elopes  occurrino",  at  a  point. 
The  data  arc  two-d       nal  and  compatible  with  the  nature  of  the  water 

i  strument  has  demonstrated  the  ability  to  function  as 
hot}:  a  laboratory  instrument  and  a  portable  field  unit.   It  is  insensi- 
tive to  wave  he]  I        ations,  thus  making  it 

>f  swell  in  either  clear 
;  crater  conditions.   The  instrument  has  an  increased  slope  measur- 
ing capacity  and  a  denonsti   ..'  ■  .   Hit;  for  daytime  field  use  not 

instrui  nts.  "  '    Errors  were  based  on  equciioia,  that 
were  not  li)  J.istic  description  of  wave- 

]  >-■   '  Y   resolute  or  :.  I     hat  of   Cox. 

For  the  tests  conducted,  the  :  ■•-. '    errors  due  to  wavelength  reso- 
lutiO!  .ted  to  be  -2.5%   at  \  =  1.33  cm,  -5%   at  X  =  -8°  Cir'> 

-II. 9^  at  X   0?  c::-  and  -20.5??  at^-  . 


CJH  AFTER  \  !  - 
:OKCLlTI)H{G  RE 


It  i?  feli  that  the  instrument  described  here  is  the  most  advanced 
wave  slope  measuring  in       t  of  its  type  to  date.  Further  improve- 
ments, such  as  greater  uj  ;  slope    ;uring  capability,  lower  tempera- 
ture-d         .  Libration  errors,  and        wavelength  resolution 
.  .,.,   •  ch.   Greater  field  ucility  may  result  from  an  added 

trv  module  and  a.  simplified  calibration  scheme  which  may 

,  ]  led  in  the  field  11  necessary.         Lruvientfs  performance  is 
expected  to  increase  markedly  by  incorporating  a  field  lens  in  place  of 

screen.   Data  pr<  'me  less  tine  consuming 

by  inter::.--,  the  lata  acquisition  system  with  a  high  speed  digital 
computer.   A" third-generation  instrument  should  provide  wave  amplitude 
informatioi  !        >n  to  the  two-dimensioi      \   slope  information, 
thus  providing  a  more  complete  time  history  of  the  surface  at  a  point. 
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